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Abstract. The profile of a critical hole in an undercooled wetting layer is deter- 
mined by the saddle-point equation of a standard interface Hamiltonian supported 
by convenient boundary conditions. It is shown that this saddle-point equation can 
be mapped onto an autonomous dynamical system in a three-dimensional phase 
space. The corresponding flux has a polynomial form and in general displays four 
fixed points, each with different stability properties. On the basis of this picture we 
derive the thermodynamic behaviour of critical holes in three different nucleation 

regimes of the phase diagram. 
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1. Introduction 



The equilibrium thickness of a wetting layer on a wall is a convenient order parame- 
ter for wetting phase transitions [1]. Fig.l shows the phase diagram for a first-order 
wetting transition in terms of temperature T and chemical potential /i [1]. Above 
the coexistence value /i c of the two-fluid bulk system the layer thickness is infinite 
whereas for h = /i — /i c < it is finite. In the limit h — > from below the layer 
thickness continuously runs to infinity above the wetting temperature T w , but has 
an infinite jump across the partial- wetting line T < T w , h — 0. A finite jump from 
a thin to a thick layer occurs when the prewetting line T p (h) is crossed from the 
region T < T p (h) to T > T p (h). This jump runs to infinity when T w is approached 
along the prewetting line and it disappears at the prewetting critical point T pc . 

A metastable wetting state can be generated by overheating a thin layer from 
T < T p (h) into a nucleation region bounded by T p (h) and an upper spinodal line 
T us (ti). The transition to the stable thick-layer configuration occurs via the random 
formation of droplets on the thin layer and growth of the supercritical droplets [2]. 
Close to the prewetting line the critical droplets have a cylindrical shape with a 
diverging radius at T p (h). This has been pointed out by Joanny and de Gennes [3] 
who have chosen the name pancake droplets for this kind of critical nuclei. 

It is also possible to undercool a thick wetting layer from T > T p (h) into a second 
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nucleation region located between T p (h) and a lower spinodal line Ti s (h). In this 
case the critical nuclei are holes in the layer which near T p {h) are mirror images of 
the pancake droplets. Close to the partial-wetting line, however, the critical holes 
have a funnel-shaped profile with a diverging depth F c but a finite inner radius 
R c at h — 0. There exists a third regime, adjacent to the wetting transition point 
T = T w , h — 0, at which F c and R c both diverge. We will refer to these regimes as 
the pre-dewetting, the partial-dewetting, and the complete-dewetting regime. 

In order to calculate the near-coexistence behaviour of F c , R c , and of the excess 
free energy E c of the critical hole in all regimes, we need some general information 
on the critical-hole profile. This is extracted from the saddle-point equation of 
a standard interface Hamiltonian [4] (which also determines the profile of critical 
droplets). As a general result this approach confirms that the size of critical nuclei 
diverges at the first-order lines h — 0, T < T w and T = T p (h), whereas it shrinks 
to zero at the spinodal lines. 

For macroscopic critical nuclei in the region h « 0, T < T w the saddle-point equa- 
tion can be mapped onto an autonomous dynamical system in a three-dimensional 
phase space. The flux of this system is polynomial and has a surprisingly rich 
fixed-point structure in the subspace of critical holes. One of these fixed points 
only exists, if the bulk dimension d of the system is smaller than some critical 
dimension d\. In this case there appear two different sets of physical trajectories 
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in the flow diagram which for T <T W correspond to critical holes at h — 0, and at 
h ^ 0, respectively. Only the latter continue to exist for d > d±, and the asymptotic 
behaviour of their trajectories determines the T, /i-dependence of F c , R c , and E c . 
For critical holes at h — the only nontrivial result of this kind is the T-dependence 
of R c . 

2. Interface model for critical holes 

Our calculations are based on the Hamiltonian [4] 



where f(x) is the local thickness of the wetting layer on the (d — l)-dimensional 
planar wall of the system. The gradient term describes long-wavelength capillary 
excitations with a stiffness constant 7. V(f) is an effective potential of the form 
shown in fig. 2, where the repulsive core simulates the wall, and V(f) — > for 
/ — > 00. The height of the minimum at / o changes with temperature, and, in 
a mean-field picture, l^(/oo) = at T = T w . Therefore, the spreading coefficient 
S = V(foo) can be used to measure the temperature distance from the wetting 
transition point [1]. 

In the full potential $(/) = V(f)—hf the minimum of V(f) at / = 00 is, for h < 0, 
shifted to a finite value /1, as in fig. 3. Along the prewetting line the two minima 
of $(/) have equal height, and they coincide at the prewetting critical point. The 
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(2.1) 



situation of fig. 3 corresponds to some point in the lower nucleation region of fig. 1. 
At the lower spinodal line fi has reached the local maximum in $(/), changing it 
into a turning point. 

Under the assumption of cylinder symmetry of critical holes, their radial profile 
f(r) obeys the saddle-point equation 

7/ "( r ) + d --l 1 f{r) = dV/df{r) - h. (2.2) 

Convenient boundary conditions for such objects are f'(r = 0) = and f(r = 
oo) = fi for h < or f'(r = oo) = for h = 0, respectively. Nontrivial solutions 
of the latter type can in fact exist in a certain range of dimensions, as discussed 
below. 

In order to prove the existence of critical holes inside the nucleation region at 
h < 0, we adopt an argument by Coleman [5]. He considered (2.2) as an equation 
of motion for a fictitious particle with position / moving in time r in a potential 
— $(/) in presence of a time-dependent friction term. At time r = the particle has 
to start with zero velocity from a position f such that asymptotically for r — > oo 
it comes to rest on top of the hill at /i, generating e.g. the critical-hole profile of 
fig.4. The existence of this marginal solution is implied by continuity from that of 
undershooting and overshooting solutions. Undershooting solutions obviously are 
produced by choosing a starting position sufficiently close to the local minimum in 
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— $(/) where the particle then eventually comes to rest. On the other hand, by 
choosing the starting position sufficiently close to the higher maximum in —$(/), 
the particle spends an arbitarily long time to reach a position /* < f\, where 
$(/*) = $(/i). This allows to neglect the friction term in the particle motion for 
/ > /* so that by energy conservation the particle will overshoot the hill at f\. 

Near interior points of the prewetting line $(/) has a double-well form which for 
the radial profile of a critical hole leads to a kink solution. In the limit T — > T p (h) 
the position R c of the turning point of the kink runs to infinity, resembling the 
behaviour of a pancake droplet. 

When some point on the line h = 0, T < T w is approached, f\ and consequently 
the critical depth F c = f 1 — f diverges. In this regime the critical-hole profile at 
macroscopic distances from the wall is determined by the asymptotic behaviour of 
V(f) for / — > oo. For long-range molecular interactions this reads 

V(f)=Af 1 -°, fora>l (2.3) 

where A is the Hamaker constant, and a = 3 or a = 4 for non-retarded or retarded 
van der Waals forces, respectively [1]. By extrapolation the macroscopic profile 
F(r) of a critical hole can then be defined as the solution of the differential equation 

7 F"(r) + d —^ 1 F'(r) = -A(a - 1)F ~ a - h (2.4) 

with the new boundary conditions F(r = R c ) = and F(r = oo) = F\ at h < or 
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F'(r = oo) = at h — 0, respectively. Undershooting and overshooting solutions 
can still be created, now by controlling the initial velocity F'(r = R c ). 

Since for r — > i? c the friction term and the field ft can be neglected in (2.4), we find 
the result 



F(r) 



fr - iL.) -+ 1 



(2.5) 



which is asymptotically valid for all ft < 0. 



For r — > oo and /i^Oa linear expansion of (2.4) in F c — F(r) leads to a Bessel-type 
differential equation. This implies the asymptotic form 



F(r) = F c 



(2.6) 



where R* = [Act (a — l)/7] _1 ^ 2 F C ( CT+1 )/ 2 , and C is an integration constant. If, as 
an approximation to the full solution of (2.4), the expressions (2.5) and (2.6) and 
their derivatives are matched at some value r = R m , it turns out that R m ~ R* 
and C is of order 1. 



For r — > oo and ft, = the left-hand side in eq. (2.4) dominates for d < d\{a) where 



dUa) = 



3(7+1 



[2.7) 



and leads to the behaviour, 



F(r) = F*D 



R, 



3-d 



(2.8) 



Here the amplitude F* = [A(a + l) 2 /8^] CT+1 Rc +1 has been adopted from the pre- 
viously derived [6] exact solution 



a + l 



of (2.4) at h = in the dimension 

d (a) = ^— - , (2.10) 
a + I 

so that D — 1 in d — d (a). Due to the boundary condition F'(r = oo) = the 
asymptotic form (2.8) implies that critical holes at h — only exist in dimensions 
d > 2. The necessity of the previously mentioned additional condition d < d\(a) 
will become clear through the following analysis. 

3. Mapping to a dynamical system 

We now define the dimensionless quantities 

rF'{r) 

X = ; — : — , 

F(r) ' 

Y = °1^±^^L_ (3 1) 

lh r 2 



2 7 F(r) ' 

and consider their dependence on the time-like variable 



t = In— (3.2) 
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where r\ is an arbitrary reference scale. Due to (2.4) we find the set of differential 
equations 

X = (3-d)X-X 2 — Y + 2Z, 

a+1 

Y = 2Y(1-^X), (3.3) 
Z = 2Z{l- l -X), 

which has the four fixed points 

-^o — Y = Z = , 

X 1 = , Y 1 = d x {(r)-d , Z 1 = 0, 

(7+1 

(3.4) 

X 2 = 3 - d , Y 2 = Z 2 = , 
X 3 = 2 , Y 3 = , Z 3 = d - 1 . 

For 1 < d < di(<r) the fixed points (3.4) are all located in the subspace X > 0, 
Y > 0, Z > 0, where the critical-hole trajectories occur (whereas F'(r) < implies 
X < for critical droplets). The subscripts of the fixed-point coordinates indicate 
the numbers of attractive principal directions of each of these points. 

In the plane Z = the fixed point Pi in fig. 5 attracts the physical trajectories 
coming from X = Y = oo which then either run to the droplet region X < or to 
the more stable fixed point P 2 . The first possibility corresponds to undershooting 
solutions of the saddle-point equation whereas the second one describes solutions 
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obeying the boundary conditions for critical holes at h — 0. In fact, the fixed- 
point value X 2 = 3 — d in connection with the definition (3.1) of X reproduces 
the asymptotic behaviour (2.8) up to an undetermined prefactor. In the limit 
d — > di(a) the fixed point Pi merges into P 2 , and in fig.5 the right section of the 
basin of attraction of P 2 collapses to zero, so that critical holes at h — no longer 
exist for d > di(a). 

If h < 0, the physical trajectories approach Pi from X = Y = Z = oo but now have 
three options to continue. Most of them either run into the droplet region or to the 
most stable fixed point P 3 , representing respectively undershooting and overshoot- 
ing saddle-point solutions where the latter behave as F(r) = —hr 2 / [2(d — 1)7] for 
r — > 00. The basins of attraction for these two sets of trajectories are separated by 
a surface which is the support of the critical-hole trajectories. 

For d < d 1 the trajectories for critical holes at h < have, contrary to our previous 
belief [7] , no chance to come close to the fixed point P 2 (which led to the erroneous 
result (9) in [7]). This is a consequence of the fact that, in agreement with (2.6), 
these trajectories for t — > 00 have to run to X = 0, Y = Z = 00. According 
to (3.3) they must, however, penetrate the plane X = 2/(er + 1) above the line 
Z = (2/(<r + 1))(Y - Yi), and for X > 2/ (a + 1) obey the condition Y < 0. This 
is incompatible with a visit of the fixed point P 2 which consequently is supposed 
to have essentially no influence on the critical- hole profile for h < 0. 
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4. Critical holes at bulk coexistence 



At bulk coexistence h = there appears an infinite set of flow lines in the X,Y- 
plane running from X = Y = oo to the fixed point P2. This means that the 
saddle-point equation (2.4) for h = has infinitely many solutions which obey the 
boundary conditions for critical holes. Only one of these solutions will, however, 
correspond to a true minimum in the variational principle 5H = 0. 

The situation can most easily be analyzed in the special dimension d = d (cr). 
Then, in terms of the variables 

v{t) = (£) ^ F(r) ' ^ ln ^' 



the saddle-point equation (2.4) assumes the form 

d 



- R2 ^ a + r4rv^ 2 

7 (a+1) 2 



(4.2) 



This again can be considered as an equation of motion for a fictitious particle, now 
without a friction term. As a consequence the particle energy 

^-^w (4.3) 

is conserved. 

Eq. (4.3) can be rewritten as 

{x-if = -—[V- X ^\y^ + 1 ( 4 - 4 ) 
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where 

X Y e 

x = y~ » y = v ' A = -— > ( 4 - 5 ) 

and e = [2/((x 2 - 1)] [(a - + l) 2 AR 2 J (A 1 )] 2,(a+1) is the maximum value of 
the potential energy in (4.3). With A taken as a parameter (4.4) analytically 
describes the full flow pattern of the system which is depicted in fig. 6. Obviously 
this pattern is symmetric with respect to the line x — 1. For x — y — 1 one obtains 
the parameter value A = 1 which consequently belongs to the separatrix running 
through the fixed point Pi (see fig. 5). This corresponds to zero kinetic energy of 
the fictitious particle when it passes the maximum of the potential in (4.3). Values 
A > 1 accordingly belong to undershooting solutions whereas for A < 1 one finds an 
infinite set of solutions obeying the critical-hole boundary conditions. The profile 
(2.9) leads to the value A = 0, i.e. a simple parabola for the corresponding flow 
line. 

For d = d (a) the energy functional H[F(r)] can be written in the form H = 
ftd-i [(^ 2 - l)AY^/(2 1 )} 2/(<T+1) u[y] where 



, > y 

u — ay 



' 1 , 1 



(4.6) 



/o *2(:r-l) 
This implies 

Su 1 _£+3 1 r 9 2 1 . 4 „. 

\V — — x z -2x--_ — -y , (4.7) 



2 c 2 

x 2 - 2x - 



a-Y 



Sx(y) A a (x-1) 2 
which shows that the variational principle Su = leads to (4.4) with A = 0. The 

special solution (2.9) therefore, in fact, corresponds to a true minimum of H . 
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5. Scaling behaviour of critical holes 

In order to calculate the quantities F c , R c , and E c for a critical hole at arbitrary- 
values of h, we use the definitions 



implied by the saddle-point equation (2.2). The second eq.(5.1) can slightly be 
simplified by use of the virial theorem which states that the potential-energy part 
in (5.1) is (3 — d)/(d — 1) times the kinetic part. This follows from the scaling 
property dH[f (ar)]/da\ a= i = which in turn is implied by the variational principle 
5H[f a (r)] = for the special set of functions f a (r) = f(ar). As a result of this 
procedure we eventually find 



where fld-i is the volume of the (d — l)-dimensional unit sphere. 

Close to the line h = 0, S < 0, we can in (5.1)-(5.3) neglect the microscopic 
increment f to F c , replace V(f) by its asymptotic form (2.3), and insert for f(r) 
the macroscopic profile F(r). This leads to the result 



$'(F C + /„) = , E c = H[f(r)} - H[F C + /„] , 



(5.1) 



and extract R c from the relation 




(5.2) 




(5.3) 



F c = 



1 



(5.4) 



(a-l)A 



13 



for h — > 0, and, in leading order, to the equations 




oo 



drr-\F'f{r) , 



(5.5) 



R c +ro 



E c 



1 




oc 



drr d - 2 {F') 2 {r) . 



(5.6) 



d-1 



Here a cut-off length r <C -R c has been introduced in order to cure the artificial 
singularity which due to the extrapolation (2.5) occurs at r = R c in the case a > 3. 

In (5.5) and (5.6) we now split off integrals running from R c + r to (1 + X)R C 
wherein the choice < A 1 allows to use (2.5). In the remaining integrals we 
transform to the scaled variables r/R* for h < and r/R c for h = 0, suggested 
by the asymptotic forms (2.6) and (2.8). This leads to a power in R* and R C: 
respectively, where the cofactors are assumed to be finite in the limit h — > and 
S — > 0. The latter assumption is supported by the analysis of chapter 3 which 
leads to expect that no further singularities will occur in addition to those implied 



On a path S = const, in the partial-dewetting regime the procedure just described 
leads for \h\ — > to a constant value of R C) and to 

E c ~ F< rP^*>W (5.7) 
with F c given by (5.4). At S — 0, i.e. in the complete-dewetting regime, we find 



by (2.5)-(2.8). 
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the behaviour 

R c ~ for a < 3 , 

R c ~ F 2 lnF c /or a = 3 , (5.8) 

R c ~ F/ -1 /or a > 1 

where again F c has the form (5.4). Moreover, in this regime we obtain 

E c ~ /or <r<3, 

£ c ~ Rt 2 \nR c for a = 3, (5.9) 
E c ~ i^ 2 /or o>3. 

For critical holes at /i = the only nontrivial result is the behaviour 

R c ~ \S\~^V for (j<3, 

i? c ~ " x In /or o- = 3 , (5.10) 
R c ~ 1 5 1 _1 /or cr > 3 

for |S| -> 0. 

In the pre-dewetting regime the asymptotic behaviour of the pancake critical holes 
will, with growing distance from T w , increasingly depend on the microscopic details 
of the potential $(/). We therefore have to go back to the relations (5.2) and (5.3), 
in which we then use the fact that f'(r) is sharply peaked at the value r = R c . 
This leads for any path S = const, to a constant value of F c at the prewetting line 
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h = h p (T), and to the relations 



R c ~ (h p {T) - h) 



-1 



(5.11) 



E c ~ (/i p (T) - ft) 



2— d 



(5.12) 



which are identical to those for pancake droplets [3]. When the wetting transition 
point is approached along the prewetting line, F c diverges as in (5.4). 

The crossover lines, separating the pre-dewetting and the partial-dewetting regime 
from the intervening complete-dewetting regime (see fig.l), are of the form \h\ ~ 
jjS'l^i. This is implied by eq. (5.5) through which the spreading coefficient S 
enters the calculation in the form S + const \h\ . 

Some of the results (5. 4), (5. 7)- (5. 12) should be accessible to experiments, although 
(5.10) strictly applies only to d < di(a). For general d the behaviour (5.10) might, 
however, be observable on a path h = const < in the partial-dewetting regime. 
The relevance of the behaviour (5.7) has been pointed out previously [8] in con- 
nection with the observed anomalously long lifetimes of undercooled wetting layers 
[9]. Direct observations of the nucleation of holes in wetting layers have also been 
reported in the recent literature [10]. 
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Figure Captions 

Figure 1: T, /z-phase diagram for a first-order wetting transition. All symbols are 
explained in the introduction. The shaded area is the complete-dewetting regime, 
separated by crossover lines from the partial- and pre-dewetting regimes. 

Figure 2: The effective interface potential V(f), where S = V(/oo). 

Figure 3: The full potential <£>(/) = V(f) — hf. Here, fx is the equilibrium thickness 
of the undercooled layer, and fi — f is the depth of the critical hole. 

Figure 4: The profile of the critical hole corresponding to the potential in fig. 4. 

Figure 5: The flow diagram of the dynamical system (3.3) with all fixed points and 
their principal directions. The shaded region is the sector in the plane X = 2/ (c+1) 
penetrated by the physical trajectories for h < 0. 

Figure 6: The flow pattern corresponding to eq. (4.4) for the limiting case a = 3 
(note that F'(r = oo) = for a = 3). 
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